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Abstract. We examine the quantization of a harmonic oscillator with inverse 
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that, for < g < 3h 2 /(8m), the system admits a U{2) family of inequivalent 
quantizations allowing for quantum tunneling through the infinite potential barrier 
at x = 0. These are a generalization of the conventional quantization applied to the 
Calogero model in which no quantum tunneling is allowed. The tunneling renders 
the classical caustics which arise under the potential anomalous at the quantum 
level, leading to the possibility of copying the profile of an arbitrary state from one 
side x > 0, say, to the other x < 0. 
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1. Introduction 



In some dynamical systems there occurs a peculiar phenomenon that classical trajec- 
tories focus on one single point after a lapse of a certain time irrespective of their initial 
conditions. This occurs typically in a harmonic oscillator, where the oscillator returns to 
the initial position periodically whatever its initial velocity is. This phenomenon underlies 
the classical caustics of geometrical optics, whose quantum version has also been studied 
earlier [1, 2]. Since the phenomenon is genuinely classical, one is tempted to ask if any 
substantial change occurs in the caustics at the quantum level. A path-integral analysis 
[3, 4] indicates, however, that for quadratic systems the focusing phenomenon remains 
essentially unchanged after quantization — it arises as a recurrence of the initial profile of 
probability distributions, accompanied by certain quantum effects [5]. 

The caustics phenomenon remains to be seen, at least classically, even when an inverse 
square potential is added to the harmonic oscillator. Indeed, the Hamiltonian, 

H M = - p i + — x * +s - 2 , (1) 

on the line — oo < x < oo admits classical solutions which exhibit periodicity for positive 
strength g > 0, implying that caustics still occur despite that the system is now non- 
quadratic. This system is in fact the quantum solvable model considered by Calogero [6] 
and, like in quadratic systems, has been argued to exhibit corresponding caustics at the 
quantum level [7]. The argument is based on the conventional quantization of the system 
(1) used for the Calogero model which assumes no probability flow to pass the singular 
(infinite) barrier at x = 0. Because of this prohibition of quantum tunneling, the system 
does not reduce to the harmonic oscillator for g — > as one naively expects. 

On the other hand, it has been known in the mathematical literature that systems 
with singularity such as the one mentioned above may have inequivalent quantizations due 
to the arbitrariness of the boundary (or connection) condition at the singularity (see, e.g., 
[8, 9] for systems on the half line, and [10, 11] for those with point interaction). In this 
paper we examine the quantization of the system (1) from this viewpoint and point out 
that, for g in the range, 

°<*<£' < 2 > 

the system admits a U(2) family of inequivalent quantizations. In contrast to the con- 
ventionl one, most of these quantizations permit quantum tunneling through the infinite 
barrier at x = 0. In particular, we find that there is a distinguished quantization possess- 
ing a smooth limit to the harmonic oscillator for g — > 0. Remarkably, in this quantization 
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the classical picture of the caustics acquires a drastic change in the quantum regime: the 
focusing occurs in two points, rather than one. We show that this anomalous phenomenon 
of quantum caustics may be used to copy the profile of an arbitrary state from one side 
x > 0, say, to the other x < 0. 

2. Classical caustics and quantum states 

Before delving into the discussion of quantization of the system (1), let us quickly 
recall how the phenomenon of the classical caustics can be observed. Let H (p, x) = E be 
the energy of the solution we are looking for. With p = mx the constant energy equation 
can readily be integrated to give the classical solution, 



The integration constant to and the energy E are altered according to the initial position 
x(ti) and velocity x(ti) chosen arbitrarily, but for any choice the particle returns to the 
original position x(U + T) = x{U) for T = kn/ui with integer k. These are the classical 
caustics appearing in the system (1). 

The quantum system corresponding to (1), too, admits exact solutions of eigenstates 
for the Schrodinger equation. Although the procedure to obtain the solutions has been 
given in various references (see, e.g., [6, 12]), we shall present here a fuller treatment paying 
a special attention to the boundary condition at x = 0. To proceed, we remove the singular 
point x = from the system to define our Hilbert space as H = L 2 (R\{0}). The boundary 
condition is then considered at the limiting points x — * ±0. For the Hamiltonian operator 
H = H(—ihd/dx,x) the Schrodinger equation for energy eigenstates reads 



For the moment we only consider the positive half line 1R+ = {x > 0}, but the negative 
half 1R_ = {x < 0} can be handled analogously using the solutions on M + . If we set 




1/2 



(3) 



Hifj n (x) 




(4) 






and choose 




(6) 
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ci = 1 + a, 

(8) 

C2 = 1 - a. 



then the Schrodinger equation (4) becomes 

+ (« + 1 " *) ^(*) " \ + 1 " A n) /„(*) = 0, A n = ^, (7) 

under the variable z = y 2 . This is the confluent hypergeo metric differential equation 
zf"(z) + (7 — z)f'(z) — af(z) = 0, whose two independent solutions are, for 7 7^ integer, 
given by f(z) = F(a, 7; z) and ^ 1_7 F(a — 7 + I, 2 — 7; z) with 7; z) being the confluent 
hypergeometric function. Thus the two independent solutions for the Schrodinger equation 
(4) are 

The general solution ip n (x) of (4) is given by a linear combination of these two, but since 
the combination may differ on the two sides K+ and H_, we have 

V>„(aO = [^V^CkD + ^V^CNJie^ + ^Vi 1 ^ (9) 

where iV^ y and AT^ are arbitrary constants and Q(x) is the Heaviside step function. 

At this point let us examine the normalizability (square integrability) of the solutions 
(8). First, since c\ > 3/2 and ci < 1/2, we observe that as x — > the solution cj)n 

approaches zero while (j)^ diverges. From J„ e dx\(j)n\x)\ 2 ~ e 2c2 for a small e, we realize 

(2) 

that </>n can be normalizable for C2 > 0. This is the case if the coupling constant g satisfies 
(2), and we confine ourselves to this case hereafter. (For the normalizability g may be non- 
positive, but for our consideration of quantum tunneling and caustics we assume g > 0.) 
Note that (2) implies 3/2 < c\ < 2 and < ci < 1/2, and this allows us to disregard the 
case 7 = integer in considering the solution of (7). Once the two independent solutions 
are admitted from the behaviour near x = 0, then the normalizability is ensured if the 
solution vanishes sufficiently fast at the infinity x = ±00. From the asymptotic behaviour 
of the confluent hypergeometric function, 

F(a, T ,z)~^e z z a -\ as \z\ - 00, (10) 

we find that the normalizability of the solutions (8) requires 

= = r(( Cl -A n )/2)r(c 2 ) 

n& n[ 2) r(( C2 -A n )/2)r( Cl y 



(ii) 
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Another condition to be imposed on the solutions is the boundary condition at the 
singular point x = 0. This is needed to ensure the continuity of the probability current 
at x = 0, which is equivalent to ensuring that the Hamiltonian H be self-adjoint. It is 
known [8, 9] that, in the presence of singularity, there can exist (at most) a U(2) family 
of self-adjoint Hamiltonians specified by corresponding boundary conditions. By means 
of the Wronskian W[ip,(p](x) = (tp(dip/dx) — {dtp / 'dx)(p)(x), which is finite even if the 
wavefunctions ip(x), <p(x) may be divergent at the singularity, the boundary conditions are 
presented as follows [13, 14] (see [15, 16] for the conditions on the line). Let ipi, ip 2 be two 
independent, real zero modes, 

Hip^x) = H<p 2 (x) = 0, W[<p u <p 2 ](x) = 1. (12) 

Given a state ip which is normalizable, we introduce the complex column vectors, 

/ w&, ^] +0 \ W[i;, H +0 \ , ^ 

Im^-oJ' \-WW,<f 2 Uj' 1 J 

defined from the boundary values W[ip, <f]±o '■= lim^^o W[ip, <f](x). The boundary con- 
dition for ip G H is then given by 

(U - I)t> + iL (U + I)t>' = 0, (14) 

where U is a U (2) matrix, / is the identity matrix, and Lq is a constant with dimension of 
length. This way a self-adjoint Hamiltonian is specified uniquely by the matrix U, which 
may hence be called the 'characteristic matrix'. 

In our case, we label n = no for which A no = in (8) and set 



V mu ( 15 ) 

<p*{*) ■= — - — 02? (M), 

so that (12) is fulfilled. Since F(a, n f; z) = 1 + O(z) as z — > 0, the boundary vectors (13) 
for the solution ifj n in (9) turn out to be 

The relations (11) and (16) then imply that the vector is proportional to and hence 
there exists a constant £ such that \P' = Thus the boundary condition (14) is now 

[(£/- J) + zL £(£/ + /)]* = (), (17) 
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and, in order to obtain a non-trivial vector we need 



det \U - I + iL £(U + I)\ = det \D - I + iL £(D + I) \ = 0, (18) 

where we have decomposed U G U (2) as U = V~ 1 DV using some SU(2) matrix V and a 
diagonal matrix D. In terms of the parameterization, 



D = ( ^_ ) , (19) 



with 9± G [0, 2tt), we find that (18) is satisfied if 



Z=-J^ or -JZ> L± = L cot^). (20) 



Substituting this back to (11), we obtain 



r((ci-A w )/2)r( C2 ) _ I h Cl -c 2 h Cl -c 2 

r((c2-A„)/2)r(ci) Vmw L+ Vmw L_ ' 1 J 

from which we determine the energy spectrum {E n = X n huj} of our system. The ratios 
N^/N^ and^ 1} /ivi 2) are determined once either L_|_ or L_ is chosen. Our result shows 
that the system permits two distinct series of eigenstates generically, one specified by L + 
and the other by L_, and this illustrates the fact that any one dimensional system which 
admits a U (2) family of self-adjoint Hamiltonians possesses a spectral family parametrized 
by two angle parameters [15, 16]. 

We shall mention a few cases where the spectrum {E n } can be obtained explicitly. 
First, if (#+,#_) = (0,0), then 1/L± = and hence (21) is fulfilled by those A n for which 
the Gamma function in the denominator has poles. This leads to E n — (2n + c 2 )/iw and the 
eigenstate given by (f>n(\x\) either on R + or 1R_ (hence each level is doubly degenerated). 
Similarly, if (6> + ,6>_) = (n, 7r), then L± = and one obtains E n = (2n + ci)Tuo and the 
eigenstate 4>n\\x\) which is also doubly degenerated. This is the case (which amounts to 
the choice U = —I) that has been considered conventionally in the treatment of the system 
(4) since the early days of Calogero [6] . 

On the other hand, if (6> + ,6>_) = (0, 7r), then 1/L+ = = L_, which means that 

(2) (2) 

there occurs two series of eigenstates, one with = ' = and the other with 

= = 0, whose eigenvalues are 

= (2n + Cl )hw, E^ = (2n + c 2 )huj, (22) 
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respectively. In particular, in the limit g — > we have c\ — > 3/2 and C2 — > 1/2, which shows 
that our system recovers the spectrum of a harmonic oscillator. A complete reduction to 
the harmonic oscillator system is realized by choosing U = o\ (where {<7j} are Pauli 
matrices), which is obtained by setting V = e l7F<T2 / 4 as well as (6> + ,6>_) = (0, n). For this 
choice, the boundary condition (17) requires N^' = —N^ \ = and hence the 
two series of eigenstates in (9) are found to be 

^\x):=N^cf>W(\x\)[e(x)-e(-x)}, 
^(x):=N^^(\x\), 

for n = 0, 1, 2..., where = [^/mw/hr(n + c s )/{(r(c s )) 2 n!}] 1 / 2 for s = 1, 2 are 

normalization constants determined so that dx \ifjn \x) | 2 = 1. The eigenfunctions 
(23) reduce exactly to those of the harmonic oscillator in the limit g — > 0, that is, ifjn^ 1 
reduces to e~ y / 2 i^2n+i(y) and ifj^ to e~ y ^ 2 H2 n (y) where H n is the Hermite polynomial 
of degree n. This in turn implies that, for other U, the system does not lead to a harmonic 
oscillator in the limit, which suggests that our system with finite g may be regarded, 
effectively, as a system that possesses a singular point interaction at x = which is hidden 
in the singularity of the potential. Under regular potentials, point interactions are known 
to admit a U(2) family of boundary conditions at the singularity, in which U = o\ provides 
the boundary condition for the 'free point', namely, no interaction there [11]. The fact that 
the smooth limit g — > to the harmonic oscillator is gained at U = o\ suggests that the 
above effective picture for the U(2) family works also for singular potentials. We also 
mention that the case U = o\ corresponds to the quantization discussed in ref.[12] which 
pointed out that the conventional quantization U = —I cannot be a perturbed harmonic 
oscillator because of the too severe physical conditions it presupposes. 

3. Quantum caustics and its anomaly 

Now that we have unconventional but perfectly admissible eigenstates arising under 
the boundary conditions specified by U, we next examine how the caustics appear at the 
quantum level. Before this, however, let us consider the possibility of quantum tunneling 
though the barrier of the potential at x = 0. In order to make our discussions clear 
and simple, we consider only the case U = a\ where the eigenstates are given by (23). To 
investigate whether or not tunneling phenomena occurs, we simply evaluate the probability 
current j(+0)(= j(— 0)) for a given arbitrary state ip. Since (23) gives our complete basis, 
we expand it as ifj(x) = ^2 n (cn ipn (x) + Cn iftn ( x ))- Then we find 

n,l 
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which shows that, since s ^ for g > 0, the probability current does flow through the 
barrier x = 0. Note that j(+0) ^ is realized for states ip consisting of both type of 
eigenstates ipn ' 1 and i>\ 2 \ and this is made possible only for g satisfying (2) and further 
for (generic) U, such as the one U = a\ we are considering, under which the two type 
of eigenstates appear. If g > 3ft 2 /8, or else if U is diagonal U = D like the conventional 
choice U = —J, we always have j(+0) = 0, disconnecting the right and left half lines, x > 
and x < 0, physically. 

Once the quantum tunneling is allowed, then the classical picture of caustics, which 
occur in the half lines independently, is no longer viable, and one is curious what in fact will 
happen quantum mechanically. To investigate this, we calculate the transition amplitude, 
the Feynman kernel K(xf,tf,Xi,ti), from the initial state of the particle staying at x = Xi 
at t = ti to the final state staying at x = Xf at t = tf. In our case (23), a straightforward 
computation (see Appendix) yields that for T := tf — ti ^ kir/uj with k = 0, 1, 2, . . . , 

Tjr . . mw ,, ..1/9 f i mu cos(uT) . 2 2 \ 

K{xf,tf;xi,ti) = . rjiA lxfXil) 1 exp -— — . [Xf+x^ 

^/ , f r f mu \xfxA \ ( mu \xtxA \1 ,^^ s 

x e(x f xi) \i a — j y± + i- a — I ; ^ U 25 

+e<- XfX .)f-i (™ fe^i! ^ +/ Mil 

where /^(z) is the modified Bessel function and a is related to (7 by (6). The last two terms 
with the factor 0(— XfXi) represent the transition allowed by the quantum tunneling. One 
can readily check that the Feynman kernel (25) reduces to that of a harmonic oscillator in 
the limit g — > 0. 

On the other hands, for T = kn/uj, we find 

K(xf,tf, Xi, ti) = (— l) fc cos(ak7r)5(xf — xi) + i(—l) h sm(akn)S(xf + xi). (26) 

The term containing 8{xf—xi) represents the quantum counterpart of the classical caustics, 
whereas the term containing S(xf + xi) represents extra caustics that arise only at the 
quantum level through the tunneling effect. We emphasize that the appearance of the 
anomalous quantum caustics is crucial to achieve the smooth reduction to the harmonic 
oscillator, since g ^ implies a — > 1/2 and hence the two terms contribute to the caustics 
of the harmonic oscillator alternately. 

In passing we note that the other limit a — > 1 is also smooth, because then the 
kernel, (25) or (26), becomes the usual one [17] (since, for (25) the last two terms with 
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Figure 1. Process of quantum copy through the caustics anomaly. At every period 
T = kn/u, a mirror image of the original profile on x > emerges on the other side 
x < 0. The relative size of the mirror image depends on a and k. 

Q(-XfXi) cancel with each other, whereas for (26) we only get S(xf —Xi)) obtained under 
the conventional quantization. This is due to the fact that, in our treatment, the second 
solution ipn \x) in (23) ceases to exist formally as a — > 1 because of the normalization 
factor N^ 2 K 

In order to see the physical consequence of the caustics anomaly, let us consider an 
initial state ip(x, U) whose density pi(x) = \ip(x, ti)\ 2 has a support only on R + . The state 
evolves according to the rule set by the kernel (25), and hence the profile will broaden and 
enter in R_ at some later time. The salient feature of the usual quantum caustics observed 
for quadratic systems is that, at tf = U + T with T = (period of caustics) x integer, the 
initial profile is reproduced completely. In our system, however, this is no longer true 
because for T = kn/uj we have the final state ip(x,tf) = f dx' K(x,tf] x' ,ti) ip(x' ,ti) with 
the density, 



This shows that, at any later periods, the profile on R + is copied as a mirror image 



image becomes exactly of the same size as the original for odd k, while for even k the 



p f (x) = \^(x,t f )\ 2 



cos 2 (a/c7r) Pi(x) + sin 2 (o/c7r) pi(—x). 



(27) 



on 1R_ (see Figure 1). In particular, when a 



3/4 (i.e., g = 5% 2 / (32m)) , the mirror 
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complete profile is reproduced on 1R_ and 1R + alternately. We note that this does not 
contradict the no-cloning theorem [18] because the two 'state spaces' on the half lines, 
L 2 (R + ) and L 2 (R_), do not comprise the entire Hilbert space by their direct product, 
H = L 2 (R \ {0}) 7^ L 2 (R + ) <g> L 2 (R_). In short, rather than making a replica of an 
arbitrary state prohibited by the no-cloning theorem, the above copying process duplicates 
a profile by the mirror image. 

Since the system discussed in this paper arises in various branches of physics, we ex- 
pect that our result will find several other applications, and to conclude we just mention 
a few. First, if one is to confine a particle with more than one channels among which the 
probability can flow like in certain nuclear states [19] or nano-devices with spin channels, 
then our quantizations may be adequate to apply. The second is the analysis of black 
holes, where our system (with and without the harmonic term) describes a particle prob- 
ing the near-horizon geometry [20, 21, 22]. Further, a straightforward extension of the 
quantizations of the n-body Calogero model (and its related solvable models) along the 
line outlined here would also enlarge the scope of the application of the model on account 
of the quantum tunneling now allowed. 

Acknowledgement: I.T. is indebted to T. Cheon and R. Sasaki for useful comments. This 
work has been supported in part by the Grant-in- Aid for Scientific Research (Nos. 10640301 
and 13135206) by the Japanese Ministry of Education, Science, Sports and Culture. 

Appendix 

In this Appendix we calculate the Feynman kernel K(x f,tf, Xi,ti) from the energy 
eigenfunctions (23). Putting T = tf — ti it is given by 

K{x f ,t r ,Xi,U) = {x f \e-^ T ' h \xi) = + S (2 \ (A.l) 

with 

oo 

S (s) = J- ^){x f ) e~^ )T s = 1, 2. (A.2) 

n=0 

To evaluate we plug (23) into (A.2) using dimensionless variables yi = yJjnujfUx^ 
yf = y^muj/hxf to find 

SW = J^(\y f yi\r- 1/2 e-^f + y^ [Q{ vm ) -Q(-y fyi )} 

. (A.3) 
x^e-^^-^L^" 1 )^)^- 1 )^), 

n=0 
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where we have used the relation between the confluent hypergeometric functions and the 
(associated) Laguerre polynomials, 



F(-n^z) = ^^L^(z). (A.4) 

Employing the standard trick T — * T(l — ie) with an infinitesimal e in (A. 3) to ensure the 
convergence of the kernel, and using the Hill-Hardy formula (see p. 189, ref.[23]), 

^ F(u + n + l) 

n=0 1 ' (A.5) 

( 1 "\ f u + v\. /2 f n (uvw) 1 / 2 \ 

= i exp -w- (uvw) V/2 I U 2^—^ , 

\l—wj \ 1—wJ \ 1 — w J 

valid for \w\ < 1, where I v (z) denotes the first kind of the modified Bessel function, we 
obtain 



SM = lim J^(| y/y< |)«i-V2 e -i(«?+«?) [Q{y yi) -Q(-y fyi )} 

X 1 _ e -i^T-e eX P I ^ \_ e-i2uT-e J (A.6) 

where we have renamed 2uTe as e for brevity. For T ^ kn/uj with fc = 0, 1, 2, . . . , we can 
take the limit e — > +0 safely to get 

'"U - ra- (ra<^ + •*>) ( , 7) 

xie^-e^,)]/.-^-^,). 

The contribution can be evaluated analogously and the result is exactly the same as 
S^ 1 ' except that c\ is now replaced by C2 and the factor [®{yfyi) — 0(— VfVi)] is removed 
in (A.6) or (A. 7). Combining the two, for T ^ kir/uj we obtain the kernel (25). 

On the other hand, for T = kn/u the kernel can be evaluated directly from (A.l). 
From the energies (A.l) and the parity ipn\— x) = (—l) s ipn\x) of the eigenstates (23), 
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we find 

oo 

s=l,2 n=0 

1 OO 

= 2 {e-^ + e-^) £ £#(*/) (V^O*))* (A.8) 

s=l,2n=0 

1 OO 

s=l,2n=0 

Using the completeness of the eigenstates and the relations c\ = 1 + a, C2 = 1 — a, we 
obtain (26). 
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